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ABSTRACT 

We employ the magnetic and velocity fields from turbulent dynamo simulations to synthesize the 
polarization of a typical photospheric line. The synthetic Stokes profiles have properties in common 
with those observed in the quiet Sun. The simulated magnetograms present a level of signal similar to 
that of the Inter-Network regions. Asymmetric Stokes V profiles with two, three and more lobes appear 
in a natural way. The intensity profiles are broadened by the magnetic fields in fair agreement with 
observational limits. Furthermore, the Hanle depolarization signals of the Sr I A4607 A line turn out 
to be within the solar values. Differences between synthetic and observed polarized spectra can also 
be found. There is a shortage of Stokes V asymmetries, that we attribute to a deficit of structuring in 
the magnetic and velocity fields from the simulations as compared to the Sun. This deficit may reflect 
the fact that the Reynolds numbers of the numerical data are still far from solar values. We consider 
the possibility that intense and tangled magnetic fields, like those in the simulations, exist in the Sun. 
This scenario has several important consequences. For example, less than 10% of the existing unsigned 
magnetic flux would be detected in present magnetograms. The existing flux would exceed by far that 
carried by active regions during the maximum of the solar cycle. Detecting these magnetic fields would 
involve improving the angular resolution, the techniques to interpret the polarization signals, and to a 
less extent, the polarimetric sensitivity. 

Subject headings: convection — line: profiles — MHD — Sun: activity — Sun: magnetic fields — Sun: 
photosphere 


1. INTRODUCTION 

In the absence of gradients of velocity and magnetic field 
in the resolution element, the polarization emerging from 
an atmosphere must be either symmetric or antisymmet¬ 
ric with respect to the central wavelength of each spectral 
line (e.g., Unno 1956; Landi Degl’Innocenti & Landolfi 
1983). The polarization observed in the solar photosphere 
does not show such symmetries. Asymmetric Stokes pro¬ 
files^ arise from the quiet Sun (e.g., Sanchez Almeida et al. 
1996; Sigwarth et al. 1999), from plages and network re¬ 
gions (e.g., Baur et al. 1981; Stenflo et al. 1984), as well 
as from sunspot penumbrae (e.g., Grigorjev & Katz 1972; 
Makita 1986; Sanchez Almeida & Lites 1992). The fact 
that such asymmetries are found even when observing at 
the limit of present resolution suggests that a rich struc¬ 
turing remains unresolved to the current observations. 

Because asymmetries carry information on the spatially 
unresolved properties of the photospheric plasma, their 
study and correct interpretation offers a chance to over¬ 
come the limitations imposed by the angular resolution. 


and to retrieve information inaccessible to direct imaging. 
Indeed, such possibility has been exploited during the last 
decade, always relying on a considerable amount of model¬ 
ing and assumptions. In broad terms, one can distinguish 
two approaches depending on the size of the unresolved 
structures that are responsible for the asymmetries. 

The first approach assumes that the unresolved photo¬ 
spheric structures are actually on the verge of being re¬ 
solved in broad-band images taken with the current in¬ 
strumentation. The smallest detectable features are of the 
order of 0.2'.'(or 150 km on the Sun); a size set by techni¬ 
cal limitations of the present solar telescopes (e.g., Bonet 
1999). Several models of this kind have been proposed to 
explain asymmetries in special cases like penumbrae, plage 
and network regions, and in the quiet Sun (e.g., Solanki 
& Montavon 1993; Grossmann-Doerth et al. 1988; Bellot 
Rubio et al. 1997; Steiner 2000). For example, Steiner 
(2000) invokes special thermal structures (temperature in¬ 
version) within magnetic regions in order to explain some 
extreme Stokes V shapes frequently observed in the quiet 
Sun. The distinctive feature of this approach is that the 


^ The term Stokes profile denotes the variation within a spectral line of any of the four Stokes parameters. Using the standard nomenclature, 
we use Stokes I to represent the intensity, Stokes Q and U for the two independent types of linear polarization, and Stokes V for the degree of 
circular polarization. Examples are given in Figure 1. 
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proposed configurations (thermal, magnetic or kinetic) are 
specific to the case under study. The fact that asymmetries 
occur everywhere is therefore difficult to explain within 
this framework. 

The second approach is more consonant with the ubiq¬ 
uity of the asymmetries in the polarization of photospheric 
lines. It assumes that the photospheric plasma is in a tur¬ 
bulent state and that small-scale structures both in veloc¬ 
ity and magnetic field are present. Order of magnitude 
estimates for the magnetic and kinetic Reynolds numbers 
in the granulation indicate that the magnetic field could 
be structured on spatial scales as small as a few kilome¬ 
ters (e.g., Schiissler 1986). This picture is consistent with 
recent progress in dynamo theory suggesting that a sub¬ 
stantial part of the magnetic field in the quiet Sun could be 
generated locally by dynamo action driven by the granu¬ 
lar flows (Meneguzzi & Pouquet 1989; Petrovay & Szakaly 
1993; Cattaneo 1999; Emonet & Cattaneo 2001). This 
second approach is also supported by certain observations 
indicating that the magnetic field is structured on spatial 
scales below the resolution limit of current telescopes (e.g., 
Berger & Title 1996; Berger et al. 1998; Sanchez Almeida 
1998). Based on these premises, it is reasonable to think 
of the asymmetries as the result of observing magnetic 
fields that vary spatially on scales much smaller than the 
mean-free-path of the photons. Interpretations of Stokes 
profiles taking into account this very fine structuring of 
the atmosphere have been carried out for several years 
(Sanchez Almeida et al. 1996). The model atmospheres 
used to fit the observed profiles consist of a collection of 
magnetic and non-magnetic components, each containing 
mild gradients to comply with the height variations of the 
mean photosphere, but interleaved in such a way as to 
produce large gradients along the line-of-sight (Sanchez 
Almeida 1997). The appealing aspect of this approach 
is that asymmetries emerge spontaneously, and indepen¬ 
dently of the details of the model; arising from correlations 
between the magnetic and the velocity fields of the various 
components (e.g., asymmetries occur if, on average, down¬ 
flowing elements have low magnetic field strength). It is 
known that atmospheres with micro-structure reproduce 
all kinds of Stokes profiles observed in the quiet Sun, in¬ 
cluding network and Inter-Network (IN) regions (Sanchez 
Almeida & Lites 2000; Socas-Navarro & Sanchez Almeida 
2002). These Micro-Structured Magnetic Atmospheres 
(MISMAs) portray a quiet Sun with large amounts of un¬ 
signed magnetic flux and very complex magnetic topology 
(very often two magnetic polarities coincide in a resolution 
element). 

In the present paper we take advantage of recent devel¬ 
opments in the numerical modeling of surface dynamos to 
understand the origin of asymmetries in the line polariza¬ 
tion. We use the magnetic and velocity fields from a set of 
numerical data generated to study the interaction between 
thermally driven turbulent convection and magnetic fields 
(Cattaneo 1999; Emonet & Cattaneo 2001). Although 
these numerical simulations were not designed for a de¬ 
tailed description of the solar photosphere, the complex¬ 
ity and ubiquity of its fields recall in many respects the 
magnetic quiet Sun inferred from the observed asymme¬ 
tries. Assuming a Milne-Eddington (ME) atmosphere for 
the thermodynamic variables we produce synthetic Stokes 


profiles. Thus, the asymmetries in the resulting profiles 
are directly related to the correlations between the veloc¬ 
ity and the magnetic field that exist in the numerical data, 
but they are decoupled from the thermodynamic variables 
of the simulation. The comparison between the synthetic 
profiles and solar data shows that the synthetic spectra are 
frequently similar to the observed ones. Such agreement 
suggests that the simulation includes some of the ingredi¬ 
ents that characterize the quiet Sun magnetic fields; in par¬ 
ticular, the correlations between magnetic field and veloc¬ 
ity at the smallest spatial scales. On the other hand, dis¬ 
crepancies between synthetic and observed Stokes profiles 
allow to identify missing ingredients and, consequently, to 
devise strategies to improve the simulations and the inver¬ 
sions. Finally, following Emonet & Cattaneo (2001), one 
can use the synthetic profiles to estimate the angular res¬ 
olution required to determine the basic properties of the 
magnetic structures present in the numerical simulation. 
This angular resolution may be of relevance to decide the 
specifications of new solar instruments. 

The work is organized as follows. The numerical data is 
briefly described in §2. The synthesis procedure is detailed 
in §3 and includes two subsections: in §3.1 we characterize 
the angular resolution of the model observations, and in 
§3.2 we calibrate the synthetic magnetograms. The main 
results of the synthesis are analyzed in §4: the variation of 
the apparent flux of the region depending on the sensitivity 
and spatial resolution of the observation (§4.1), the asym¬ 
metries of the Stokes V profiles (§4.2), and the magnetic 
broadening of the intensity profiles (§4.3). Hanle depolar¬ 
ization signals to be expected for Sr I A4607 A are worked 
out in §5. §6 discusses the diameter of the ideal telescope 
needed to spatially resolve the simulations. Finally, the 
implications of the present work are discussed in §7. 

2. DESCRIPTION OF THE SET OF MHD DATA 

The numerical data set used in the present work is part 
of a series of numerical simulations to study the generation 
and interaction of magnetic fields with turbulent convec¬ 
tion. The simulations are based on a idealized model de¬ 
scribing a layer of incompressible (Boussinesq) fluid with 
constant kinematic viscosity p, thermal diffusivity k and 
magnetic diffusivity rj. The boundary conditions are pe¬ 
riodic in the horizontal directions, and impenetrable and 
stress free with constant temperature and zero horizontal 
magnetic field along the upper and lower boundaries of 
the computational domain. As is common in this kind of 
simulations, the unit of length is the vertical extent of the 
layer, the unit of time is the thermal diffusion time across 
that layer and the magnetic intensity is expressed as an 
equivalent Alfven speed. With these units the particular 
dynamo solutions used below is defined by the following di¬ 
mensionless parameters: aspect-ratio of the computational 
domain 10 x 10 x 1, Rayleigh number 5 x 10^ and Prandtl 
numbers vjn = 1, p/r/ = 5. The numerical resolution 
was 512 X 512 x 97 collocation points. In this regime, the 
solution corresponds to a state of vigorous convective tur¬ 
bulence characterized by a kinetic and magnetic Reynolds 
numbers of Re = 200 Re^ — 1000, respectively. The 
resulting flow is strongly chaotic and acts as an efficient 
dynamo generating an intermittent magnetic field with no 
mean flux and a magnetic energy roughly 20% of the ki- 
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netic energy. For the present study, we use one time step of 
the dynamo evolution that is well into the statistically sta¬ 
tionary regime. At this epoch, the rms velocity of the fluid 
is u « 200 with a corresponding turnover time^ of 1/100. 
A more complete description of the simulation procedure, 
and further details about the numerical solutions can be 
found in Cattaneo (1999), Emonet & Cattaneo (2001) and 
Cattaneo et al. (2001). 

3. SPECTRAL SYNTHESIS 

The numerical simulations described above were not de¬ 
signed specifically for spectral synthesis, rather to study 
dynamo action. In order to achieve the high magnetic 
Reynolds numbers (Re^ = 1000) necessary for dynamo ac¬ 
tion to manifest itself, and given finite computer resources 
(lO’^-lO® grid points), we found it necessary to resort to 
the Boussinesq approximation. The latter is valid provided 
that the vertical extent of the layer is much smaller than 
the pressure, density, and temperature scale heights (see 
e.g. Spiegel & Veronis 1960). As a consequence, the ther¬ 
modynamic variables from the numerical solutions have 
modest variations across the layer, and cannot be used to 
calculate the opacity and source functions needed for the 
synthesis. Consequently, we address the synthesis prob¬ 
lem ignoring the thermodynamic variables from the simu¬ 
lations, but keeping the variations in the velocity (u) and 
magnetic field (b). Although u and b are then decoupled 
from the density and temperature distributions, the result¬ 
ing synthetic profiles are useful to study the properties of 
the emerging line polarization that depends, primarily, on 
the magnetic and velocity structure. In the remainder of 
this section we explain the process followed to generate the 
synthetic profiles. 

The first step is to translate the dimensionless velocity 
and magnetic data of the simulations into units that are 
practical for the synthesis. The spatial scale is fixed by as¬ 
suming that the typical size of one granule is about 1000 
km. In the statistically stationary regime the computa¬ 
tional domain is about 5 granules wide. Thus the hori¬ 
zontal size of the domain corresponds to 5000 km, with a 
corresponding grid resolution of approximately 10 km in 
the horizontal, and 5 km in the vertical. Magnetic field 
strengths B and velocities U readily follow from the en¬ 
ergy equipartition arguments. B and U are proportional 
to the dimensionless magnetic strength b and velocity u 
through the factors /„ and ft, 

B = fbb, U = fuu. (1) 

Equipartition between kinetic energy density and mag¬ 
netic energy density is reached in dimensionless units when 

b% = ( 2 ) 

or equivalently when 

B%/{^^) = pUlj2 , (3) 

p being the density. The previous equation, together with 
equations (1) and (2), lead to 


Thus, if one associates the rms speed u ~ 200 with the 
rms velocity fluctuations observed in the solar granulation 
{7 = 2 km s“^ (e.g., Durrant 1982; Topka & Title 1991), 
then 

fu = 10“^ km s~^ , (5) 

and for a typical density at the base of the photosphere (3 
xl0“^ g cm“^; Maltby et al. 1986), equations (4) and (5) 
render 

/b - 2 G . (6) 

Except for a few exploratory calculations described in §4.2 
and §5, the scaling factors in equations (5) and (6) are used 
for the rest of the work. Notice that reasonable deviations 
from this scaling do no affect our results in a substantial 
way (a point addressed in §4). 

Eollowing our previous assumption, the absorption and 
emission terms of the radiative transfer equations (e.g., 
Beckers 1969; Wittmann 1974; Landi Degl’Innocenti 1976) 
are not based on the thermodynamic variables of the nu¬ 
merical data set. Instead, we resort to a Milne-Eddington 
(ME) synthesis which assumes a constant line absorp¬ 
tion and a linear source function (e.g., Unno 1956; Landi 
Degl’Innocenti 1992). The ME assumption is routinely 
used for magnetic field diagnostics (e.g., Socas-Navarro 
2001) since it yields an analytic expression for the emerg¬ 
ing polarized spectrum that only depends on a few param¬ 
eters. The analytic solution is used to fit observed spectra 
and retrieve information from them (e.g., Skumanich & 
Lites 1987). The ME synthesis hides all the unknown ther¬ 
modynamic properties of the atmosphere in five parame¬ 
ters, namely, the line absorption coefficient, two parame¬ 
ters that define the (linear) source function, the damping 
coefficient, and the Doppler width. In our syntheses, we 
adopt the values characteristic of the Fe I line at A6302.5 
A in the quiet Sun (Sanchez Almeida et al. 1996, Table 
1, including a Doppler width of 40 mA or 1.9 km s“^). 
This magnetically sensitive spectral line is often used for 
magnetic studies, including those of the Advanced Stokes 
Polarimeter (Elmore et al. 1992). Working with the same 
line allows us to compare our results directly with its ob¬ 
servations. The use of observed ME thermodynamic pa¬ 
rameters to represent the synthetic line provides realism 
to the part of the synthesis that we do not obtain directly 
from the numerical simulation. The emerging spectra will 
have the observed equivalent widths, line widths, etc^. 

The standard ME assumption considers an atmosphere 
having a uniform magnetic field (e.g., Unno 1956; Landi 
Degl’Innocenti 1992). Such assumption is clearly at odds 
with our numerical data in which the magnetic diffusive 
scale is of the order of 500 x Rem km ~ 15 km, i.e., 
three grid points in the vertical direction. The field is 
therefore not constant over the range of heights where the 
typical photospheric lines are formed (say 100 km). A real¬ 
istic spectral line synthesis would have to take into account 
the contributions from many layers in the simulation. The 
brute-force approach would be to assign an optical depth 
to each horizontal plane in the simulation and then syn¬ 
thesizing the spectrum by direct numerical integration of 


fb = /„(47rp)i/2 ^ (4) 

^ We define the turnover time to be twice the vertical crossing time based on the rms speed u. 

^ Two independent results of the synthesis reflect this realism. First, the constant to calibrate magnetograms worked out in §3.2 agrees with 
those derived from real solar spectra. Second, the mean equivalent width of our synthetic intensity profiles is 71 mA, i.e., close to the value 
observed in the quiet Sun (some 80 mA; e.g., Moore et al. 1966). 
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the radiative transfer equations for polarized light. Here, 
instead, we resort to a different strategy that a) is simpler 
and much faster, b) reduces the number of free parameters 
of the synthesis (no need to specify the optical depths of 
all the different planes of the atmosphere) and, c) produces 
spectra very close to those obtained by the brute-force ap¬ 
proach. We employ a multi-component ME atmosphere, 
as defined by Sanchez Almeida et al. (1996; §2.1). We 
assume that not all individual planes of the simulation 
contribute to the emerging spectrum, and consider only a 
sub-set whose thickness is of the order of the photon mean- 
free-path. We take into account that the planes are opti¬ 
cally thin so that the emerging spectrum depends only on 
the average along the line-of-sight of the emission and the 
absorption (the so-called MISMA approximation; Sanchez 
Almeida et al. 1996). As argued below, these assumptions 
are not unreasonable, and reduce the radiative transfer 
calculation to deciding how many planes contribute to the 
emerging spectrum. Subsequently, the average absorption 
matrix is computed, and used jointly with the analytic 
solution for the regular ME synthesis (Sanchez Almeida 
et al. 1996). 

In our case we choose the 21 uppermost planes of the 
simulation, excluding the upper boundary. According to 
the scaling described in the previous paragraph, this cor¬ 
responds to approximately 100 km in the atmosphere, a 
typical range for the formation of a photospheric line. We 
further assume that the numerical simulation is placed at 
disk center, i.e., with the line-of-sight along the vertical 
direction. 

Since the individual planes of the simulation are indeed 
optically thin (5 km), using the average emission and ab¬ 
sorption instead of solving the full radiative transfer should 
be a good approximation. We checked this assumption by 
comparing spectra computed by direct integration of the 
radiative transfer equation with those resulting from our 
approximation. Figure 1 contains one example. The full 
synthesis is performed assuming the 21 different planes of 
the ME synthesis to be equi-spaced in and spanning from 
logr = 0 to logr = — 1. (The symbol t stands for the con¬ 
tinuum optical depth.) This layer is periodically repeated 
to complete the atmosphere from logr = 10 to logr = —3. 
The synthesis of 1000 randomly chosen spectra shows rela¬ 
tive differences between the two kinds of syntheses of only 
a few per cent. This deviation is negligible for the analy¬ 
sis that we carry out in the present work. The choice of 
number of planes to be used for the spectral synthesis is 
a more delicate matter. Using 26 instead of 21 planes we 
find differences between spectra of up to 25%. These dif¬ 
ferences are due to the strong variability of the magnetic 
conditions, an uncertainty that also affects any other way 
of synthesizing the spectra. For example, in the case of 
the brute-force approach the signals depend on the optical 
depths arbitrarily assigned to each plane of the numerical 
data. 

3.1. Spatial resolution: seeing and telescope diameter 

Studying the consequences of observing the simulations 
with a finite spatial resolution is one of the key objectives 
of the present work. We use the term seeing to denote all 
the effects that may reduce the spatial resolution (from 
genuine seeing to optical aberrations of the instruments). 


Here the effects of seeing are modeled by smearing the 
2D maps of the Stokes profiles with Gaussian functions. 
Note that the 2D convolution is carried out at each indi¬ 
vidual wavelength of each Stokes parameter. The amount 
of smearing is characterized by the FWHM (Full Width 
Half Maximum) of the Gaussian function. We take this 
particular point spread function for the sake of simplicity, 
however, it accurately describes long-integration-time at¬ 
mospheric seeing (e.g, Roddier 1981, §4.5). In one of the 
sections (§ 6 ), we measure seeing in terms of the diameter, 
D, of an ideal telescope that provides a given angular reso¬ 
lution. Considering that the Airy disk of an ideal telescope 
has a FWHM ~ A/D, then 

D = 13 cm (A/6302) (FWHM/725)-\ (7) 

where FWHM is expressed in km and A in A. We have 
used a scale of 725 km/", corresponding to our target at 
1 AU. 

3.2. Longitudinal Magnetograms 

Longitudinal magnetograms are images showing the de¬ 
gree of circular polarization in the flank of a spectral line 
(e.g., Babcock 1953; Landi Degl’Innocenti 1992; Keller 
et al. 1994). They can be produced using relatively simple 
and compact instruments and therefore have been widely 
used in solar magnetic field studies. A significant frac¬ 
tion of what we know about quiet Sun magnetic fields has 
been derived from them (e.g., Wang et al. 1995, and ref¬ 
erences therein). Longitudinal magnetograms are usually 
calibrated in units of magnetic field strength (G). Since 
we wish to compare the polarization signals with observed 
magnetograms, calibrated magnetograms must be pre¬ 
pared out of the synthetic Stokes profiles. We follow a pro¬ 
cedure that mimics the observational process. First, the 
circular polarization signals are averaged in wavelength, 
this is because magnetograms are often obtained through 
a rather broad color filter. We employ a running box filter, 
100 mA wide, centered in the blue wing of the line at 50 
mA from the line center. (These values are comparable to 
those employed in real measurements; e.g., Yi & Engvold 
1993; Berger & Title 2001.) Then the wavelength aver¬ 
aged circular polarization signal Va is calibrated to yield 
the so-called longitudinal magnetic flux density Dios 7 i-e., 

Dios = Goal 14. (8) 

The calibration constant Goal is evaluated using the mag¬ 
netograph equation and the properties of the line in which 
the magnetograph operates; explicitly, 

Go-J = (9) 

Here, the symbol Aq stands for the central wavelength of 
the line, 5 e// is the effective Lande factor, and la cor¬ 
responds to the intensity profile smeared with the same 
color filter used to produce V/. The constant c equals 
4.67 X 10 A ^ G For the atomic parameters of 
Fe I A6302.5 A, and using the mean Stokes I profile over 
a complete snapshot, we find Gcai — 5850 G when Va is 
given in units of the continuum intensity. This synthetic 
calibration constant is in good agreement with the values 
that can be found in the literature for this line (e.g., Yi & 
Engvold 1993; Lites et al. 1994; Socas-Navarro & Sanchez 
Almeida 2002 ). 
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With several additional simplifying assumptions^, it can 
be shown that 



where the integrals extend to the resolution element S, and 
Bz is the component of the magnetic field vector pointing 
towards the observer. This identity provides the rationale 
to denote magnetogram signals as longitudinal magnetic 
flux densities. (According to equation [10], Bios represents 
the magnetic flux per unit area.) 

4. RESULTS 

In this section we analyze the synthesis of the 512 x 512 
Stokes spectra calculated from the dynamo simulation de¬ 
scribed in §2. We only use the numerical solution at one 
instant since the purpose of the work is to study the kind 
of spectra characterizing the stationary state of the simu¬ 
lation (studies of the time dependent behavior is deferred 
for later). We place the simulation at the solar disk center 
so that the vertical direction follows the line-of-sight. 

Figure 2 shows the magnetogram that results when ap¬ 
plying equation (8) to the synthetic Stokes V spectra 
(right panel). The effect of observing with a 0'.'5 seeing 
- representative of the best angular resolution achieved 
at present (e.g., Berger & Title 2001), is also shown on 
the left panel. Ignoring the underlying substructure, an 
observer would identify a number of magnetic concentra¬ 
tions in the smeared magnetogram (e.g., points o, b and c). 
However, these are difficult to associate with single struc¬ 
tures once the underlying substructure is acknowledged. 
Rather, the points with enhanced signal in the 0"5 mag¬ 
netogram represent locations where the convective flows 
are continuously advecting magnetized plasma to balance 
the plasma that constantly disappears in the downdrafts 
(see also § 3.2 in Emonet & Cattaneo 2001). The fact 
that the residual polarization signal shows up suggests the 
presence of a large-scale structure of the advecting velocity. 
Indeed, detailed analysis of the numerical data reveals the 
existence of a mesogranular flow (Cattaneo et al. 2001). 

Except for a few test calculations in §4.2, the spectra 
discussed in this section correspond to a scaling of the di¬ 
mensionless magnetic field and velocity given by equations 
(5) and (6). However, syntheses using /„ = 0.02 km s“^ 
and /{, = 2.5 G were also tried. We found that the degree 
of circular polarization increases linearly as ft, whereas 
a larger /„ enhances the line asymmetries. On the other 
hand, the kind of general trends and properties that we 
analyze here do not depend on the precise value of the 
scale factors. 

4.1. Flux density versus angular resolution 

The amount of magnetic flux and energy in the numer¬ 
ical data is far larger than that detected in the Sun as 
IN fields. This difference can be understood as the result 
of two important factors that hinder the detection of the 
weak signals emerging from the simulation. First, the mag¬ 
netic fields are highly disorganized so that the polarization 
signals tend to cancel as the angular resolution deteriorate 
(Emonet & Cattaneo 2001). Second, most synthetic sig¬ 
nals are extremely low, i.e. at and below the sensitivity 

^ Weak magnetic field that is constant along the line-of-sight, Stokes 
DeglTnnocenti & Landi DeglTnnocenti (1973); Jefferies & Mickey (1991). 


of the present instrumentation. Consequently a large frac¬ 
tion could not be detected at present. In this section we 
use our synthetic lines to explore these effects. We find 
that once realistic angular resolution and sensitivity are 
taken into account, the simulations are in good agreement 
with the observations. 

Figure 3 shows the mean unsigned signal in the magne¬ 
togram of Figure 2 as a function of the seeing (i.e., mean 
I Bios I over the snapshot versus FWHM of seeing). We con¬ 
sider various sensitivities of the magnetograph: moderate 
(20 G), good (5 G) and very good (0.5 G). We account for 
the limited sensitivity by setting to zero all those points in 
the smeared magnetogram where the signal is below the 
hypothetical observational threshold. Figures 3b is iden¬ 
tical to Figure 3a except that it has been normalized to 
the mean longitudinal field in the simulation (i.e., mean 
[B^l ~ 51 G, where the average considers all the points of 
the simulation that we use). According to the standard in¬ 
terpretation (equation [10]), this is the parameter that one 
retrieves from a magnetogram. The normalization helps 
visualizing the fraction of real signal that remains in the 
magnetogram for a given angular resolution and sensitiv¬ 
ity. There are several features in these two figures that 
deserve comment. 

Even with perfect sensitivity and maximum angular res¬ 
olution, one detects only 80 % of the existing flux. The 
cancellation is mostly due to the averaging along the line- 
of-sight caused by the radiative transfer. At maximum 
angular resolution, the flux in the magnetogram depends 
little on the sensitivity since most of the signals exceed 
20 G. 

The decrease in signal strength as the angular resolution 
deteriorate is severe. For the typical 1" angular resolution, 
the detectable signals are only 10 % of the original ones 
(Fig. 3b). This estimate is optimistic since it holds true 
if the sensitivity is good; should the latter be moderate, 
only traces of the original signals remain (1% for 20 G 
sensitivity). 

Figure 3a includes values for solar IN magnetic flux den¬ 
sities observed by various authors. The level of detected 
flux density agrees well with the predictions of the simula¬ 
tions once the angular resolution has been taken into ac¬ 
count. Note that the observed flux densities are more than 
a factor of ten smaller than the signals in the original mag¬ 
netogram. We should not overemphasize the agreement, 
since the observational points are rather uncertain (they 
come from inhomogeneous sources, with different sensitiv¬ 
ities and based on disparate techniques: see below). How¬ 
ever, two conclusions can be drawn. First, intense yet tan¬ 
gled magnetic fields like those in the numerical simulation 
are compatible with the present solar observations. Sec¬ 
ond, if our spatially fully resolved synthetic spectra were 
emitted by the Sun, they would produce the degree of 
polarization that we detect on the Sun with the present 
instrumentation. 

The observations presented in Figure 3a required trans¬ 
formation from the numbers quoted in the original works 
to the quantities plotted in the figure. Here we briefly 
describe these transformations for the sake of complete¬ 
ness. (The symbols accompanying the citations corre- 

I independent of the magnetic field, and others; see, e.g., Landi 
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spend to those used in Fig. 3a.) Socas-Navarro & Sanchez 
Almeida (2002; bullet sign •) deduce 10 G mean flux den¬ 
sity for the quiet Sun fields, however, they mention that 
the apparent flux density decreases by a factor 2.4 if it 
is estimated from magnetograms. This renders 4.2 G for 
a 1" angular resolution, a figure derived from the cut off 
frequency in the Fourier domain of the continuum inten¬ 
sity image (Sanchez Almeida & Lites 2000). The same 
procedure is used to estimate the angular resolution of 
the spectra in Gollados (2001; square □). The mean flux 
density, 3.4 G, has been directly provided by the author. 
Wang et al. (1995; asterisk *) point out, explicitly, a 1.65 
G flux density for their angular resolution of 2". Stolpe 
& Kneer (2000; times sign x) find some 3.5 G mean flux 
density (private communication), to which we associate 
an angular resolution of two pixels or 1'.'4. Lin & Rim- 
mele (1999; triangle A) do not directly give a mean flux 
density for their measurements. However, their Figure 3 
contains the magnetic fluxes of individual measurements, 
and the authors provide the scale factor between the mag¬ 
netic flux and the magnetic flux density (5x 10^^ Mx G“^). 
With this conversion, the mean flux density of these ob¬ 
servations is about 4 G. The magnetic features occupy 68 
% of the surface, since the rest remains below the sensi¬ 
tivity threshold. Gonsequently, the mean flux density over 
the whole surface results 0.68 x 4 G. Grossmann-Doerth 
et al. (1996, Fig. 1; plus sign -f) provide the distribu¬ 
tion of Stokes F Fe I A5250 A signals found in the quiet 
Sun, being the average signal about 2 x 10“^ (in units of 
the continuum intensity). The authors also point out a 
calibration constant from Stokes V signal to flux density. 
It gives 2 G for the mean observed signal. Since these 
signals fill 40% of the solar surface, the mean flux den¬ 
sity is about 0.8 G. We assign to these data an angular 
resolution twice the sampling interval (some 2"3). Finally, 
Dominguez Gerdeha et al. (2002; diamond o) have recently 
found 17 G mean flux density in a 0”5 angular resolution 
magnetogram taken with Fe I A6302.5 A. 

Figure 4 provides a different illustration of the similari¬ 
ties between the synthetic magnetogram and the real Sun. 
It shows a real magnetogram of the quiet Sun obtained 
in the blue wing of Fe I A6302.5 A (Sanchez Almeida 
& Lites 2000, Fig. 1), i.e., using the observational setup 
that we have tried to reproduce (see §3.2). A patch of 
the real magnetogram with the size of the numerical sim¬ 
ulation has been replaced with the synthetic one (Fig. 4 
left). Except for the location of the inset, which we place 
in an IN region, we have not used free parameters to pro¬ 
duce the combined magnetogram. The angular resolution 
has been chosen to match the observation (1"), whereas 
the degree of polarization comes directly from the synthe¬ 
sis. Despite the absence of fine-tuning, it turns out to be 
extremely difficult to distinguish the magnetogram with 
the artificial inset (Fig. 4 left) from the real one (Fig. 4 
right). In other words, the synthetic magnetogram con¬ 
tains structure with a spatial distribution similar to the 
observed one, and with a degree of polarization which fits 
within the observed range. 

4.2. Asymmetries of the Stokes V profiles 

In the absence of gradients of magnetic and velocity 
fields within the resolution element, the Stokes profiles 


have to obey well defined symmetries (see §1). This kind 
of symmetries are never found in the quiet Sun, instead 
profiles frequently show extreme asymmetries (Sanchez 
Almeida et al. 1996; Grossmann-Doerth et al. 1996; Sig- 
warth et al. 1999). Because the magnetic field in the 
numerical simulation has a highly intermittent structure 
(Gattaneo 1999; Emonet & Gattaneo 2001), we expected 
the resulting synthetic Stokes profiles to be asymmetric. 
In the following, we compare our synthetic profiles with 
observed ones and study in a statistical sense their differ¬ 
ences and similarities. We proceed by first classifying the 
types of Stokes V profiles produced by the simulation. The 
2.6 X 10® Stokes V profiles in the snapshot were assorted 
using a cluster analysis algorithm identical to that em¬ 
ployed with real data by Sanchez Almeida & Lites (2000; 
§3.2). The procedure identifies and groups profiles with 
similar shape, irrespective of their degree of polarization 
and polarity, since it employs Stokes V profiles normalized 
to the largest blue peak (see Sanchez Almeida & Lites 2000 
for details). 

The classification is summarized in Figure 5. For each 
different class, the averaged Stokes V profile over the en¬ 
semble is plotted. The mean is taken after each individ¬ 
ual profile has been multiplied by the sign of its largest 
blue peak. This avoids cancellation of different polarities. 
The mean profiles are numbered (upper left corner in each 
panel) according to the percentage of profiles in the simu¬ 
lation that are elements of their class, with ffO being the 
most probable and #17 the least probable. The proba¬ 
bility to find a profile of a given class is indicated in the 
upper right corner of each panel. Note the various possi¬ 
bilities. From virtually anti-symmetric Stokes V (implying 
no asymmetry), to profiles with one (e.g., #11 and #13) 
or three lobes (e.g., #12 and #14). All these asymmetries 
are produced by gradients of magnetic and velocity fields 
along the line-of-sight. The peak polarization is consider¬ 
able, some 1% in units of the continuum intensity. This 
has to be the case to yield the large flux density shown in 
Figure 3a for perfect angular resolution. Another detail 
worthwhile noting is the balance between the number of 
profiles having a large blue lobe and those whose principal 
lobe is the red one. 

The line shapes in Figure 5 are difficult to compare with 
observed profiles since observations of the quiet Sun have 
much lower angular resolution. Figure 6 shows another 
classification having all the features explained above, ex¬ 
cept that the simulation data have been smeared with a 1" 
seeing, typical of the real observations. First, the polar¬ 
ization signals are reduced by one order of magnitude with 
respect to the original profiles; compare them with those in 
Figure 5. Second, new more complicated line shapes have 
arisen as the result of the large spatial smearing (radiative 
transfer smoothes over only 100 km, whereas the spatial 
smearing does it over 725 x 725 km^). Stokes profiles hav¬ 
ing all these very asymmetric shapes are indeed observed 
in the real quiet Sun (see Sanchez Almeida et al. 1996, 
Fig. 12; Sigwarth et al. 1999, Fig. 7; Sanchez Almeida 
& Lites 2000, Fig. 4, 5 and 6; Sigwarth 2001, Fig. 4). 
This qualitative agreement is again a notable feature of 
the simulation, since there was no obvious a priori reason 
to expect it. 

However, despite such general qualitative agreement, 
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one can find quantitative differences between the synthetic 
and the observed profiles. As it happens with the fully re¬ 
solved profiles (Fig. 5), the number of synthetic profiles 
having a principal blue lobe and those with a principal 
red lobe is similar (e.g., the pairs #12 and #18; #31 and 
#34; #33 and #35). Such balance is not present among 
the observed profiles, where the blue lobe usually domi¬ 
nates (e.g., Sanchez Almeida & Lites 2000, Fig. 4; Sig- 
warth et al. 1999, Fig. 12). Another qualitative difference 
with observations concerns the profiles that are most fre¬ 
quently obtained in the simulation. They show almost 
no asymmetry (see classes #0 to #7). Contrariwise, the 
observational counterpart have a well defined asymmetry 
characterized by a large blue lobe (similar to those ob¬ 
served in plage and enhanced network regions). This lack 
of significant asymmetry can be traced back to the origi¬ 
nal syntheses (Fig. 5, classes #0 to #4), and therefore to 
the variations along the line-of-sight of the magnetic field 
and velocity. The probable cause is discussed in the next 
paragraph. 

Figure 7 summarizes the kind of variations along the 
line-of-sight existing in the numerical simulation. Several 
definitions are required before it can be interpreted. The 
variations are described using standard statistical param¬ 
eters, namely, the line-of-sight mean value /^, the line-of- 
sight standard deviation Sf^j, and the line-of-sight corre¬ 
lation coefficient 

riz 

f ij — ^ ' fijk, 

fc=l 


<5/y= (/yfc-/y)2 


1/2 


( 11 ) 


= [ifijk - fij)igi3k - 9i,)] j P7y ■ 


The arrays fijk and gijk may represent any component of 
the fields, and their indexes vary according to the position 
in the horizontal plane {i and j) and in the vertical di¬ 
rection (fc). The symbol stands for the number of grid 
points in the vertical direction. We should bear in mind 
that the line-of-sight averages defined in equations (11) 
depend on the horizontal coordinates. We need to char¬ 
acterize the typical properties of these line-of-sight mean 
values for each range of mean longitudinal magnetic field 
strength. For this purpose we define the average (h) and 
the dispersion Ah of the quantity h among all those points 
in the simulation with a given line-of-sight mean longitu¬ 
dinal magnetic field # i.e., 

nx 
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Note that the symbol h may represent any of the statistical 
parameters in equation (11) (/, 5f or (fg)), and nx and 
Hy correspond to the grid points in the two horizontal di¬ 
rections. The bin size e has to be chosen to guarantee hav¬ 
ing enough points per bin. Following these definitions, we 
can now interpret Figure 7. The different quantities rep¬ 
resented are plotted as function of the line-of-sight mean 
magnetic field. 

The first panel in Figure 7 shows the variation of (Uz) 
and (SUz). The limited dispersion of the vertical velocities 
provides an explanation for the predominance of profiles 
with small asymmetries. The dispersion must be of the 
order of the line width to produce a substantial modifi¬ 
cation of the line shape. The typical standard deviation 
of the vertical velocities along the line-of-sight turns out 
to be between 0.2 and 0.3 km s“^ {{SUz), Fig. 7a, the 
solid line), whereas the line widths are of the order of 2 
km s“^ (see §3). The magnetic field itself is probably not 
responsible for the moderate asymmetries since its varia¬ 
tions along the line-of-sight are large. For example, Figure 
7b shows the standard deviation of the mean longitudinal 
magnetic field (SBz), which is frequently larger than the 
absolute value of Bz- In fact, the variations of the longitu¬ 
dinal magnetic fields are so important that Bz very often 
changes sign along the line-of-sight (Fig. 7c). Note that 
none of these statements on large magnetic field gradients 
apply to the intrinsically strong fields, a case discussed 
separately in the next paragraph. Above, say, 200 G, the 
variations of field strength are mild and the longitudinal 
field maintains a constant sign along the line-of-sight (Figs. 
7b and 7c). Concerning the balance between the number 
of asymmetries towards the blue and towards the red in 
Figures 5 and 6, it is probably due to the lack of a def¬ 
inite sign for the correlation between magnetic field and 
velocity. Works on the asymmetries in plage and network 
regions repeatedly indicate the need for a negative corre¬ 
lation to account for the observed preponderance of the 
Stokes V blue lobe, explicitly. 


{Uz\Bz\)<0, (16) 

(see, e.g., Solanki & Pahlke 1988; Sanchez Almeida et al. 
1988; Sanchez Almeida et al. 1989; Grossmann-Doerth 
et al. 1988; Grossmann-Doerth et al. 1989; Sanchez 
Almeida 1998). Such condition is satisfied when downflows 
and magnetic fields are spatially separated, i.e., when the 
strongest downflows tend to occur in weakly magnetized 
plasma. (Note that Uz > 0 corresponds to downflows.) 
One can see in Figure 7d that the correlation between \Bz\ 
and Uz has a well defined negative value only for the largest 
field strengths. Since theses points represent a small frac¬ 
tion of the synthetic profiles, they contribute very little to 
the classification in Figure 5 which, consequently, shows 
no obvious preference for a blue asymmetry. 

Those patches in the simulation with the largest field 
strength show a clear negative correlation fulfilling the cri¬ 
terion in equation (16) (see Fig. 7d). Do they produce the 
observed Stokes V profiles with a main blue lobe? They 
do not, since the asymmetry of the profiles emerging from 
concentrations of intense field strength are minimal. This 
fact can be understood using Figures 7a and 7b, which re¬ 
veal gradients of both magnetic field and velocity too small 
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for the requirements described in the previous paragraph. 
Despite this apparent disagreement, the key ingredients to 
yield the right shapes are already present in the simula¬ 
tion. If one artificially increases the gradients of magnetic 
field and velocity already existing in the simulation, large 
asymmetries similar to the observed ones automatically 
show up. Figure 8 contains synthetic profiles emerging 
from the three more intense magnetic concentrations in the 
snapshot (labeled as a, b and c in Fig. 2, left). Note that 
they already have the Stokes V asymmetries that charac¬ 
terizes network and IN regions. (Figure 8d includes one of 
these observed Stokes V profiles for reference, namely, one 
network profile in Sanchez Almeida & Lites 2000). In or¬ 
der to produce theses new synthetic profiles with enhanced 
asymmetries, we increased the intensity of the flow with 
respect to the scaling in §3 by using /„ = 0.04. In addi¬ 
tion, we increased the variations of magnetic field strength 
by averaging all the absorption and emission over Oi'5 of 
the simulation (i.e., over all the points in a box of this 
size). The qualitative agreement between synthetic and 
real profiles is remarkable. 

The above considerations lead to an interesting ques¬ 
tion. How should the numerical simulations be modified 
in order to produce asymmetries closer to the observed 
ones? First, the dispersion of velocities at the smallest 
scales must be increased, which implies a decrease in the 
kinematic viscosity. Second, magnetic and non-magnetic 
regions should be even more intermittent to strengthen 
the correlation (16). This can be achieved by decreasing 
the magnetic diffusivity, which would both increase the 
tangling of magnetic field lines and the dispersion of field 
strengths existing in the intense concentrations. 

4.3. Broadening of the Intensity profiles 

One of the observational constraints on the existence of 
a complex and tangled magnetic field in the solar photo¬ 
sphere comes from the work by Stenflo & Lindegren (1977) 
(see also Unno 1959; Stenflo 1982). If a tangled magnetic 
field exist, it has to broaden the spectral lines of the solar 
spectrum according to their magnetic sensitivities (i.e., ac¬ 
cording to their effective Lande factors). All other things 
being the same, those with larger sensitivity should be 
slightly broader. Stenflo & Lindegren (1977) looked for 
such effect in the solar unpolarized spectrum with no suc¬ 
cess. From the error budget of the measurement the au¬ 
thors set an upper limit to the field strength of the existing 
fields (Stenflo & Lindegren 1977, equation [12]), 

Sapp < 140 G. (17) 

In this section we analyze if the magnetic fields in the sim¬ 
ulations produce line broadenings compatible with such 
observational upper limit. This exercise allows to address 
the question of whether magnetic fields as intense as those 
in the simulation may exist on the Sun, and still remain be¬ 
low the present observational detection limit. We already 
know that the degree of circular polarization of the sim¬ 
ulation stays well within the observational bounds (§4.1). 
Here we address the question from a different perspective, 
using a totally different observational constrain that de¬ 
pends on the magnetic fields in a intrinsically different 
way. 

Our synthetic Stokes I spectra have an excess of broad¬ 
ening caused by the presence of magnetic fields. Figure 


9 shows the difference between the mean Stokes I profile 
produced by the region, and the mean profile produced 
when the syntheses are repeated with no magnetic field, 
but keeping everything else identical. The magnetic pro¬ 
file is broader and shallower at the line center, producing 
a residual with three lobes. Is this extra broadening com¬ 
patible with the observational limit (17)? A detailed mod¬ 
eling of the procedure employed by Stenflo & Lindegren 
(1977) is clearly beyond our possibilities, since it requires 
the synthesis of hundreds of spectral lines with different 
temperature and magnetic field sensitivities. Fortunately, 
the essence of the procedure is simple. If two lines are 
identical except for their magnetic sensitivity, the small 
excess of width Aw can be directly related to an apparent 
magnetic field strength Happ according to the rule, 

Sapp = Ja/ (Aw/w). (18) 

The scale factor <5 depends of the wavelength, the differ¬ 
ence of Lande factor, and the mean line width w. For two 
lines with the wavelength and strength of Fe I A6302.5 A, 
one magnetic and another one non-magnetic, the scaling 
factor turns out to be 

5-2.1x10^0, (19) 

which follows from the equations (7) and (8) in Stenflo & 
Lindegren (1977). The estimate of the apparent magnetic 
field strength using equations (18) and (19) is a matter 
of determining the excess of broadening Aw/w associated 
with the residuals in Figure 9. Following Stenflo & Linde¬ 
gren (1977), the two mean intensity profiles where fitted 
using Gaussian functions, i.e., 

/c — / (X exp [ — (A/w)^], (20) 

where A stands for the wavelength relative to the line cen¬ 
ter and Ic represents the continuum intensity. Then the 
difference between the widths w for the syntheses with and 
without magnetic field directly yields 

Aw/w - 1.1 X 10-2, (21) 

or, using equations (18) and (19), 

Sapp - 220 G. (22) 

Note that the Gaussian fits reproduce fairly well the dif¬ 
ference between the two synthetic profiles (see the dashed 
line in Fig. 9). 

The field strength of our synthetic profiles (equation 
[22]) apparently come into conflict with the observational 
limit in equation (17). Should the inconsistency be real, 
it points out an excess of magnetic fields in the numerical 
simulations as compared to the solar case (an excess of 
magnetic field strength, area coverage of the fields, etc.). 
However, the marginal discrepancy is probably not signif¬ 
icant in view of the uncertainties affecting both the syn¬ 
theses and the observational limit. Actually, the similarity 
between the observational limit and the predicted width 
should be understood as real chance to test whether com¬ 
plex tangled fields like those in the numerical simulations 
are present in the solar photosphere. A slight refinement 
of the currently available diagnostic tools (e.g., Stenflo & 
Lindegren’s technique) should be able unambiguously to 
confirm them or discard them. We return to this point in 

§7. 
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5. HANLE SIGNALS PRODUCED BY THE DYNAMO 
SIMULATIONS 


Up to now we only consider polarization signals gen¬ 
erated by Zeeman effect. However, solar turbulent mag¬ 
netic fields have been inferred using Hanle effect signals®. 
They indicate the presence of turbulent magnetic fields 
in the upper photosphere-lower chromosphere with a field 
strength between 5 and 60 G (e.g., Faurobert-Scholl et al. 
1995; Bianda et al. 1999). These strengths may seem 
too low compared to those assumed in this work, with a 
mean value larger than 100 G (e.g., §4.3). Therefore, we 
felt compelled to estimate the Hanle signals expected from 
the dynamo simulation and compare them with observed 
values. Once more we face the question of whether mag¬ 
netic fields similar to those in the simulations may exist 
and still produce observable effects within solar values. 

A complete Hanle effect synthesis similar to that carried 
out for the Zeeman effect is clearly beyond the scope of this 
paper. Fortunately, one can estimate the level of Hanle sig¬ 
nals for one of the typical lines used in Hanle effect based 
diagnostics. Considering the Sr i line at A4607 A, the de¬ 
polarization produced by a magnetic field can be expressed 
as 






i + tA 


1 + 47 ^. 


(23) 


where QjQo is the ratio between the observed linear po¬ 
larization Q and the polarization expected if there were 
no magnetic field Qo- The symbol 7 // parameterizes the 
magnetic field strength of the turbulent field B, 


JH = B/Bh- 


(24) 


The normalization factor Bh scales linearly with the ra¬ 
diative transition rate plus the depolarizing collision rate. 
The relationship (23) is an approximation which works 
well for this particular line in standard quiet Sun model 
atmospheres (see Faurobert et al. 2001, where one can also 
find the dependence of Bh on the density and temperature 
of the atmosphere). For additional details see, e.g., Landi 
Degl’Innocenti (1985) and Trujillo Bueno & Manso Sainz 
(1999). 

In case that the magnetic field strength is not unique 
but has a distribution of values, the mean signal Q/Qo 
turns out to be equal to the mean depolarization factor 
considering the distribution of the field strengths Wb (see 
Landi Degl’Innocenti 1985, §3). Note that Wb no longer 
varies with B as B ^ Bh {'Jh ^ 1), therefore the aver¬ 
age depolarization Wb bears no information on the large 
field strengths that may exist in the distribution. In other 
words, the mean signal will always be biased towards weak 
field strengths (e.g., Trujillo Bueno 2001, Fig. 3). This 
fact could reconcile the large fields in the dynamo simu¬ 
lations with the observed Hanle depolarization signals. In 
order to check such possibility, we evaluate the mean Wb 
for the snapshot studied in the paper. The transition and 
collision rates required to compute Bh are determined, ac¬ 
cording to the prescription in Faurobert et al.(2001, §3), 
using the quiet Sun densities and temperatures given by 
Maltby et al. (1986). We only consider the photospheric 
layers. The expected signals, including all the magnetic 


fields in the snapshot describing the distribution of B, are 
represented in Figure 10. The depolarization changes with 
height due to the variation of temperature and density in 
the model atmosphere, which modify Bh- Figure 10 in¬ 
cludes the range of Hanle signals observed close to the disk 
center, namely, Q/Qo = 0.50±0.18 (quoted by Faurobert- 
Scholl et al. 1995, from the measurements of Stenflo 1982 
at an heliocentric angle whose cosine is /x = 0.81). The ob¬ 
servational signal has been chosen to be as close as possible 
to the disk center to sample deep photospheric layers. The 
measurement of depolarization is ascribed by Faurobert- 
Scholl et al. (1995) to a range of heights between 200 
km and 400 km. Although these layers are still too high 
for the range of heights that we assign to the simulation 
in the Zeeman syntheses (the base of the photosphere), 
the depolarization turns out to be within the bounds set 
by observations (see Fig. 10, the solid line). The agree¬ 
ment improves taking into account that the turbulent field 
strengths are expected to decrease with height in the at¬ 
mosphere (e.g. Faurobert-Scholl et al. 1995). The dashed 
line in Figure 10 has been computed assuming the field 
strengths to be a factor of two smaller than those in the 
original snapshot. In short, our simplified estimate in¬ 
dicates no obvious inconsistency between the simulations 
and the observed Hanle effect depolarization signals for 
Sr I A4607 A. 

6. WHAT SPATIAL RESOLUTION IS NECESSARY TO 
RESOLVE THE NUMERICAL DATA? 

Answering this question may be relevant to design of 
the next generation solar telescopes (Rosner & Beckers, 
2001, private communication). For the sake of simplicity, 
we adopt the magnetic flux as the physical parameter to 
be determined. It is routinely measured using standard 
instrumentation, and it suffers a severe bias due to insuf¬ 
ficient spatial resolution (§4.1). The fraction of magnetic 
flux in the simulation that one still detects in synthetic 
magnetograms is used to quantify the required angular res¬ 
olution. This information is contained in Figure 3b, which 
we plot in a slightly different way in Figure 11. Abscis¬ 
sae and ordinates have been interchanged, and the angular 
resolution is also presented in terms of the diameter of the 
ideal telescope achieving the angular resolution (see equa¬ 
tion [7]).This representation suits our present purposes. 
The curve is shown for two different sensitivities, i.e., when 
(virtually) all polarization signals are detected, and when 
only those larger than 20 G are above the noise level. The 
following conclusions follow from Figure 11: 

1. Irrespective of the resolution, one only detects 
80 % of the existing flux. The residual 20 % 
cancellation is due to the smearing along the 
line-of-sight produced by the radiative transfer. 
Surpassing this upper limit is not a question of 
increasing the telescope size but rather it demands 
improving the diagnostic techniques used to 
measure the magnetic flux. 

2. Detecting 50 % of the flux would require an angular 
resolution of 110 km or a telescope of 85 cm. 


® The Hanle effect is a purely non-LTE phenomenon whose details and subtleties are still a subject of active research. Roughly speaking, the 
magnetic field modifies the polarization of the light that we detect after scattering (see, e.g., Landi DeglTnnocenti 1992; Stenflo 1994). 
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Detecting another 20 % additional magnetic flux 
demands resolving 40 km and so a telescope some 
three times larger (2.3 m). 

3. Points #1 and #2 refer to observations without 
noise. If one considers a detection threshold of 
20 G (which corresponds to a degree of polarization 
of some 0.3 %, according to the calibration in 
§3.2), then only 70 % of the flux present in the 
simulation can be detected. On the other hand 
detecting 50 % implies an angular resolution of 
70 km or a diameter of ~ 1.3 m. 

4. The slope of the solid line in Figure 11 drastically 
changes when trying to detect more than, say, 

65 % of the original flux. Going beyond this point 
requires a large increase of telescope diameter for 
a limited increase of additional signal. This 65 % 
may represent an optimal compromise between 
resolution and telescope size, and it corresponds to 
60 km on the Sun or a diameter of 1.5 m. 

7. DISCUSSION 

The dynamo simulation by Gattaneo (1999) (see also 
Emonet & Gattaneo 2001) produce magnetic fields whose 
structure resembles in many ways the quiet Sun magnetic 
fields (§1). Although the simulations were not designed for 
a realistic description of the solar conditions, the quan¬ 
titative comparison with the quiet Sun that we under¬ 
take is useful for a variety of reasons. It allows to judge 
whether, and to what extent, the turbulent dynamo pro¬ 
vides a paradigm to describe the quiet Sun magnetism. 
It allows to identify, study and understand the difficulties 
and biases faced by the current observational techniques 
when applied to very complex fields. It helps identifying 
physical ingredients that are missing in the simulations, an 
exercise helpful to guide future numerical work. Finally, it 
may suggest the type of technical developments needed to 
measure the properties of a magnetic field with the com¬ 
plexity present in the numerical data. Keeping in mind 
all these reasons, we synthesize the polarization emerging 
from the simulation. Specifically, we choose a snapshot 
of the time series representative of the stationary regime, 
which renders some 2.6 x 10® individual spectra. The 
assumptions and limitations of our Milne-Eddington ap¬ 
proach to the synthesis are discussed in §3. They let us 
calculate the polarization produced by a magnetically sen¬ 
sitive spectral line similar to Fe I A6302.5 A. The synthetic 
spectra have been analyzed in the light of observations of 
quiet Sun magnetic fields. The comparison reveals simi¬ 
larities and differences. In addition, it provides some hints 
and caveats to keep in mind when interpreting observa¬ 
tions. These three aspects of the synthetic spectra are 
discussed next. 

The magnetograms produced by the numerical data are 
in good agreement with those observed in the Sun (§4.1, 
Figs. 3a and 4). The agreement, however is only achieved 
after a severe canceling of signals produced by the poor 
spatial resolution of the current observations. More than 
90 % of the unsigned magnetic flux existing in the numeri¬ 
cal data does not appear in magnetograms with l"angular 
resolution (Fig. 3b). In other words, it is this 2 to 10 
% residual that agrees with the observed signals. The 


Stokes V profiles emerging from the simulation are fre¬ 
quently very asymmetric, in qualitative agreement with 
observations (see references in §1). In particular, Stokes V 
profiles with three and more lobes result from the existence 
of two opposite polarities in the resolution element. Part 
of these asymmetries are produced by gradients of mag¬ 
netic held and velocity along the line-of-sight (Fig. 5), 
but the Stokes V profiles also owe much of their shapes 
to variations across the line-of-sight within the resolution 
element (cf. Figs. 5 and 6). We studied in §4.3 the ex¬ 
cess of broadening of the synthetic intensity profiles due 
to the presence of a magnetic field. The additional broad¬ 
ening is found to be close to, but typically in excess of, the 
upper limit of 140 G set by Stenflo & Lindegren (1977). 
The uncertainties of the measurement could easily explain 
this discrepancy, therefore, we understand this marginal 
disagreement as an invitation to revisit the work by Sten¬ 
flo & Lindegren (1977) and improve the sensitivity by a 
factor two. This should be enough clearly to confirm or 
discard a solar magnetic field with the features present in 
the numerical data. So far only Zeeman signals of typical 
photospheric lines have been mentioned. We also estimate 
in §5 the Hanle depolarization to be expected if the simu¬ 
lation is placed at various heights in the photosphere. All 
the uncertainties notwithstanding, the depolarization sig¬ 
nals for Sr I A4607 A turn out to be within the observed 
bounds. 

A quantitative analysis of the Stokes V asymmetries re¬ 
veals real discrepancies between the synthetic spectra and 
the quiet Sun. Most of the synthetic profiles show mild 
asymmetries, well below the mean observed values. More¬ 
over, the synthetic spectra do not contain the clear ob¬ 
served tendency for the Stokes V blue lobe to dominate. 
Such trend is present even in the largest polarization sig¬ 
nals, which trace big concentrations of magnetic held. We 
believe that the cause of the discrepancy is twofold. First, 
the dispersion of velocities between nearby pixels is too 
small. Second, the spatial separation between strongly 
magnetized and unmagnetized plasmas is too large. These 
two factors minimize the asymmetries in the intergranular 
lanes, despite the fact that the kind of correlation between 
magnetic field and velocity already existing in the simu¬ 
lation produces the right asymmetry. Arguments for such 
explanation are given in §4.2. As a support of this view, 
we artificially increase the dispersion of velocities and mag¬ 
netic fields existing in the simulation to synthesize spectra 
in three particularly large magnetic concentrations. The 
resulting Stokes V profiles look very much like those ob¬ 
served in the quiet Sun network (synthetic and observed 
profiles are represented in Fig. 8). If we have correctly 
identified the origin of the differences, what must be mod¬ 
ified in the numerical simulation in order to produce asym¬ 
metries closer to the observed ones? One would need to 
increase the dispersion of velocities and magnetic fields 
at the smallest spatial scales, keeping the relationship be¬ 
tween them already existing in the simulation. This could 
be attained by increasing the Reynolds numbers of the 
simulations. 

There is another discrepancy between the synthetic and 
the observed profiles which has not been mentioned yet. 
The analysis of Stokes V shapes observed in the quiet Sun 
reveals that kG held strengths are common (see Sanchez 
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Almeida & Lites 2000, Socas-Navarro & Sanchez Almeida 
2002). There are also weaker fields (e.g., Lin & Rimmele 
1999; Collados 2001), but the frequency of kC is certainly 
larger than that present in the simulation. This differ¬ 
ence can be readily pin down to the incompressibility of 
the simulation (§2), which hampers the strong evacuation 
of the plasma, an ingredient always associated with the 
existence of kG fields in the photosphere. 

We have found several similarities between the synthetic 
spectra emerging from the dynamo simulation and obser¬ 
vations of quiet Sun magnetic fields. Such (sometimes un¬ 
expected) agreement seems to point out that the simula¬ 
tion already capture some of the essential features charac¬ 
terizing the quiet Sun magnetic fields. These arguments 
enables us to proceed by analogy and propose the exis¬ 
tence of magnetic fields on the Sun similar to those in the 
numerical simulations. Several interesting consequences 
follow from this assumption: 

1. The amount of flux inferred by conventional 
magnetograms depends very much on the angular 
resolution and the sensitivity of the measurements. 
Good or very good sensitivity is mandatory to 
detect any fluxes with the canonical l"angular 
resolution. For example, the structures in the 
numerical simulation would be very difficult to 
detect in magnetograms like those provided by the 
MDI instrument on board the SOHO spacecraft, 
with a sensitivity of some 15 G (e.g., Liu & 
Norton 2001). According to Figure 3b, only 1% 
of the original magnetic flux is detectable at 
this sensitivity level. On the other hand, the 
sensitivity is not so critical upon improvement 

of the angular resolution. In the limit of perfect 
angular resolution the signals become of the order 
of 50 G, corresponding to a degree of circular 
polarization of 1% (see Fig. 5). 

2. The amount of unsigned flux existing in the 
simulation is also very large in absolute terms. 
Assuming that the full solar surface were covered 
by magnetic fields like those in the numerical 
simulations, the associated total unsigned flux 
would be of the order of 3.1 x 10^^ Mx. This 
figure is a factor 4.5 times larger than the total 
flux detected at solar maximum with conventional 
techniques (Rabin et al. 1991, Fig. 5). Should 
such large amounts of hidden flux exist in the 
photosphere, it would have to interact with 

the other classical manifestations of the solar 
magnetism (active regions, solar cycle, coronal 
heating, etc). This point in particular deserves 
further investigations. 

3. We discuss in §6 the telescope size required 
to detect the magnetic structures observed in 
the simulations. This is a preliminary estimate 
which certainly has to be revised and refined. 


(For example, we do not address the question of 
whether the 60 km optimum spatial resolution is 
set by the photon mean-free-path or the diffusion 
processes included in the MHD simulation.) 
However one conclusion of this preliminary study 
stands out. A sizeable fraction of the magnetic 
flux is not detected due to the radiative transfer 
smearing along the line-of-sight. Improving such 
bias is not a matter of increasing the telescope 
size and polarimetric sensitivity. It only depends 
on interpreting the observed polarization using 
techniques that account for the line-of-sight 
smearing. Keeping in mind the complexity of the 
fields (e.g., the longitudinal magnetic field almost 
always changes sign along the line-of-sight; see Fig. 
7c), the use of inversion techniques that assume 
optically thin fluctuations of the magnetic field 
seems to be a reliable selection (e.g., the MISMA 
inversion code in Sanchez Almeida 1997). 

4. The Stokes asymmetries contain information on the 
magnetic structuring at spatial scales smaller than 
the resolution element, an information unique and 
valuable (see §1). The syntheses show that both 
variations along the line-of-sight and across the 
line-of-sight are responsible for the asymmetries. 
Gonsequently, the inversion techniques aimed 
at understanding the origin of the asymmetries, 
and therefore at measuring the properties of the 
photospheric magnetic fields, have to incorporate 
the two aspects of the spatial variations. 

The importance of the above conclusions depends to some 
extent, on the degree of realism of the simulations. How¬ 
ever, the tendencies that they imply are probably correct; 
e.g., the existence of large amounts of solar magnetic flux 
missing in conventional measurements, or the need for non¬ 
standard inversion techniques to determine the properties 
of the quiet Sun magnetic fields. These problems are im¬ 
portant enough to deserve a close follow-up. 

The interest to synthesize the spectra of turbulent dy¬ 
namo simulators arose from conversations with T. Bogdan 
back in 1998. Both R. Rosner and J. Beckers pointed out 
the importance of estimating the telescope size required 
to spatially resolve the simulation. Thanks are due to F. 
Kneer, F. Stolpe, and M. Gollados for providing the flux 
densities cited in §4.1, and to an anonymous referee for 
helpful comments. Thanks are also due to J. Stenflo for 
clarifications on the Hanle effect based measurements of 
the turbulent magnetic field strength. Discussions with 
R. Manso and J. Trujillo Bueno were extremely helpful 
to elaborate §5. This work has been partly funded by 
the Spanish Ministry of Science and Technology, under 
projects AYA2001-1649 and DOES 95-0028-G. Two of the 
authors FG and TE were partially supported by NASA 
grant NAG5-10831 at the university of Ghicago. 



12 


Sanchez Almeida, Emonet & Cattaneo 


REFERENCES 


Babcock, H. W. 1953, ApJ, 118, 387 

Baur, T. G., Elmore, D. E., Lee, R. H., Querfeld, C. W., Sz Rogers, 
S. R. 1981, Sol. Phys., 70, 395 
Beckers, J. M. 1969, Sol. Phys., 9, 372 

Bellot Rubio, L. R., Ruiz Cobo, B., Sz Collados, M. 1997, ApJ, 478, 
L45 

Berger, T. E., Lofdahl, M. G., Shine, R. A., Sz Title, A. M. 1998, 
ApJ, 506, 439 

Berger, T. E., k Title, A. M. 1996, ApJ, 463, 365 
Berger, T. E., k Title, A. M. 2001, ApJ, 553, 449 
Bianda, M., Stenflo, J. O., k Solanki, S. K. 1999, A&A, 350, 1060 
Bonet, J. A. 1999, in ASSL, Vol. 239, Motions in the Solar 
Atmosphere, ed. A. Hanslmeier k M. Messerotti (Dordrecht: 
Kluwer), 1 

Cattaneo, F. 1999, ApJ, 515, L39 

Cattaneo, F., Lenz, D., k Weiss, N. . 2001, ApJ, 563, L91 
Collados, M. 2001, in ASP Conf. Ser., Vol. 236, Advanced Solar 
Polarimetry - Theory, Observations, and Instrumentation, ed. 

M. Sigwarth (San Francisco: ASP), 255 
Dominguez Cerdeha, I., Kneer, F., Sanchez Almeida, J. 2002, 
ApJ, submitted 

Durrant, C. J. 1982, in Landolt-Bornstein, Vol. 2a, Astronomy • 
Astrophysics k Space Research, ed. K. Schaifers k H. H. Voigt 
(Berlin: Springer-Verlag), 96 
Elmore, D. F., et al. 1992, Proc SPIE, 1746, 22 
Emonet, T., k Cattaneo, F. 2001, ApJ, 560, L197 
Faurobert, M., Arnaud, J., Vigneau, J., k Frish, H. 2001, A&A, 
378, 627 

Faurobert-Scholl, M., Feautrier, N., Machefert, F., Petrovay, K., k 
Spielfiedel, A. 1995, AkA, 298, 289 
Grigorjev, V. M., k Katz, J. M. 1972, Sol. Phys., 22, 119 
Grossmann-Doerth, U., Keller, C. U., k Schiissler, M. 1996, A&A, 
315, 610 

Grossmann-Doerth, U., Schiissler, M., k Solanki, S. K. 1988, A&:A, 
206, L37 

Grossmann-Doerth, U., Schiissler, M., k Solanki, S. K. 1989, A&:A, 
221, 338 

Jefferies, J. T., k Mickey, D. L. 1991, ApJ, 372, 694 
Keller, C. U., Deubner, F.-L., Egger, U., Fleck, B., k Povel, H. P. 
1994, AkA, 286, 626 

Landi Degl’Innocenti, E. 1976, A&AS, 25, 379 
Landi Deghlnnocenti, E. 1985, Sol. Phys., 102, 1 
Landi Degl’Innocenti, E. 1992, in Solar Observations: Techniques 
and Interpretation, ed. F. Sanchez, M. Collados, k M. Vazquez 
(Cambridge: Cambridge University Press), 71 
Landi Degl’Innocenti, E., k Landi Degl’Innocenti, M. 1973, 

Sol. Phys., 31, 299 

Landi Degl’Innocenti, E., k Landolfi, M. 1983, Sol. Phys., 87, 221 
Lin, H., k Rimmele, T. 1999, ApJ, 514, 448 

Lites, B. W., Martinez Pillet, V., k Skumanich, A. 1994, Sol. Phys., 
155, 1 

Liu, Y., k Norton, A. A. 2001, Mdi measurement errors: the 
magnetic perspective. Technical Report SOI 01-144, Center for 
Space Science and Astrophysics, Stanford University 
Makita, M. 1986, Sol. Phys., 106, 269 

Maltby, P., Avrett, E. H., Carlsson, M., Kjeldseth-Moe, O., Kurucz, 
R. L., k Loeser, R. 1986, ApJ, 306, 284 


Meneguzzi, M., k Pouquet, A. 1989, J. Fluid Mech., 205, 297 
Moore, C. E., Minnaert, M. G. J., k Houtgast, J. 1966, The Solar 
Spectrum from 2935 A to 8770 A , NBS Mono. 61 (Washington: 
NBS) 

Petrovay, K., k Szakaly, G. 1993, AkA, 274, 543 
Rabin, D. M., Devore, C. R., Sheeley, N. R., Harvey, K. L., k 
Hoeksema, J. T. 1991, Solar interior and atmosphere, ed. A. N. 
Cox, W. C. Livingston, k M. S. Matthews (The University of 
Arizona Press), 781 

Roddier, F. 1981, in Progress in Optics, ed. E. Wolf, Vol. 19 
(Amsterdam: North-Holland Publishing Co.), 281 
Sanchez Almeida, J. 1997, ApJ, 491, 993 
Sanchez Almeida, J. 1998, ApJ, 497, 967 

Sanchez Almeida, J., Collados, M., k del Toro Iniesta, J. C. 1988, 
AkA, 201, L37 

Sanchez Almeida, J., Collados, M., k del Toro Iniesta, J. C. 1989, 
AkA, 222, 311 

Sanchez Almeida, J., Landi Degl’Innocenti, E., Martinez Pillet, V., 
k Lites, B. W. 1996, ApJ, 466, 537 
Sanchez Almeida, J., k Lites, B. W. 1992, ApJ, 398, 359 
Sanchez Almeida, J., k Lites, B. W. 2000, ApJ, 532, 1215 
Schiissler, M. 1986, in Small Scale Magnetic Flux Concentrations 
in the Solar Photosphere, ed. W. Deinzer, M. Knolker, H. H. 
Voigt (Gottingen: Vandenhoeck k Ruprecht), 103 
Sigwarth, M. 2001, ApJ, 563, 1031 

Sigwarth, M., Balasubramaniam, K. S., Knolker, M., k Schmidt, 
W. 1999, AkA, 349, 941 

Skumanich, A., k Lites, B. W. 1987, ApJ, 322, 473 
Socas-Navarro, H. 2001, in ASP Conf. Ser., Vol. 236, Advanced 
Solar Polarimetry — Theory, Observations, and Instrumentation, 
ed. M. Sigwarth (San Francisco: ASP), 487 
Socas-Navarro, H., k Sanchez Almeida, j. 2002, ApJ, 565, 1323 
Solanki, S. K., k Montavon, C. A. P. 1993, A&:A, 275, 283 
Solanki, S. K., k Pahlke, K. D. 1988, AkA, 201, 143 
Spiegel, E. A., k Veronis, G. 1960, ApJ, 131, 442 
Steiner, O. 2000, Sol. Phys., 196, 245 
Stenflo, J. O. 1982, Sol. Phys., 80, 209 

Stenflo, J. O. 1994, Solar Magnetic Fields, ASSL 189 (Dordrecht: 
Kluwer) 

Stenflo, J. O., Harvey, J. W., Brault, J. W., k Solanki, S. K. 1984, 
AkA, 131, 333 

Stenflo, J. O., k Lindegren, L. 1977, A&:A, 59, 367 
Stolpe, F., k Kneer, F. 2000, A&A, 353, 1094 
Topka, K. P., k Title, A. M. 1991, Solar interior and atmosphere, 
ed. A. N. Cox, W. C. Livingston, k M. S. Matthews (The 
University of Arizona Press), 727 
Trujillo Bueno, J. 2001, in ASP Conf. Ser., Vol. 236, Advanced 
Solar Polarimetry - Theory, Observation, and Instrumentation, 
ed. M. Sigwarth (San Francisco: ASP), 161 
Trujillo Bueno, J., k Manso Sainz, R. 1999, ApJ, 516, 436 
Unno, W. 1956, PASJ, 8, 108 
Unno, W. 1959, ApJ, 129, 375 

Wang, J., Wang, H., Tang, F., Lee, J. W., k Zirin, H. 1995, 

Sol. Phys., 160, 277 
Wittmann, A. 1974, Sol. Phys., 35, 11 
Yi, Z., k Engvold, O. 1993, Sol. Phys., 144, 1 



Polarization of lines from turbulent dynamo simulations 


13 



- 400-200 0 200 400 - 400-200 0 200 400 

Wavelength [mA] Wavelength [mA] 

Fig. 1. — Approximation used to synthesize the polarized spectrum. The solid lines show Stokes 7, Q, U and V profiles from a point of the 
simulation synthesized according to the MISMA approximation employed along the text. The dashed lines represent a full integration of the 
radiative transfer equations. Typical relative deviations are of the order of a few per cent. Wavelengths are in mA and the Stokes profiles 
have been normalized to the continuum intensity. 



Fig. 2. — Magnetogram of the snapshot of the numerical simulation analyzed in the work. It shows both the original one at full resolution 
(right), as well as the result of smearing the magnetogram with a 0(^5 seeing (left). The white gauge indicates 1" on the Sun. The labels a, b 
and c point out three entities that may be identified as single magnetic concentrations in the 0(^5 seeing magnetogram. The scale of grays of 
the two magnetograms is independent since the seeing diluted signals (left) would remain almost unnoticed if scaled to those in the original 
magnetogram (right). 
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Seeing [Mm] Seeing [Mm] 

Fig. 3. — Mean flux density as a function of the spatial resolution of an hypothetical observation. The three different curves represent 
different sensitivities of the observation (0.5 G, the solid lines; 5 G, the dotted lines; 20 G, the dashed lines), (a) The flux density in the 
synthetic magnetogram is given in G. The symbols represent real observations of IN fields. For the equivalence between symbol and reference, 
see §4.1. (b) The flux density has been normalized to the true mean flux density, i.e., the mean value of \Bz\ in the simulation. It helps 
understanding the large fraction of missing flux for the sensitivity and angular resolution of a typical observation (say, 5 G and l'^). Note the 
20 % reduction for perfect angular resolution, which is produced by the radiative transfer smearing along the line-of-sight. The arrow on the 
abscissae points out 1". 
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Fig. 4.— Synthetic magnetogram embedded in the real magnetogram of a quiet Sun region: inset within a box in the left image. The right 
image shows the original magnetogram, including the box for reference. The bar in the bottom corresponds to 25 000 km on the Sun, i.e., the 
linear dimension of a typical network cell. Note how the synthetic magnetogram fits in smoothly within the real magnetogram. It produces 
the right polarization for IN fields and it also has the proper spatial scales. 
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Fig. 5.— Types of Stokes V profiles produced by the numerical simulation. All plots share a common scale which is indicated at the bottom 
of the figure. The label of each profile includes a number for cross-reference, and the percentage of profiles in the simulation belonging to 
the category. The profiles have been taken from the original simulation, with no spatial smearing. The large plus sign of each plot indicates 
the origin of the horizontal and vertical scales, shown at the bottom of the figure. The Stokes V profiles are normalized to the continuum 
intensity. 
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|Bj[G] |BJ[G] 

Fig. 7. — Variations along the line-of-sight of the longitudinal magnetic field and the longitudinal velocity. We represent mean values among 
all those points in the simulation having the line-of-sight mean longitudinal magnetic field given in the abscissae {Bz \ see equations [11] to [15] 
for definitions), (a) Mean vertical velocity (the dashed line) as well as the mean standard deviation among the velocity fluctuations along the 
line-of-sight (the solid line). {Uz > 0 corresponds to downfiows or redshifts.) (b) Mean standard deviation of the magnetic field fiuctuations 
along the line-of-sight. The oblique straight line y = x \s included for reference, (c) Percentage of points where the longitudinal magnetic 
field changes sign along the line-of-sight. (d) Correlation coefficient for the fluctuations of magnetic field and velocity along the line-of-sight. 
The dashed line corresponds to the mean value, {{Uz\Bz\))^ whereas the shaded region shows one standard deviation about this mean value, 
{{Uz\Bz\)) di A(C/ 2 |-B 2 |). Velocities and magnetic fields are in km s~^ and G, respectively. The correlation coefficients are dimensionless. 
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Fig. 8. — (a), (b), and (c) correspond to the synthetic Stokes V profiles resulting form the average of the patches a, b and c in Figure 2. 
(d) Observed profile shown for reference. It is the class 3 profile in Figure 4 of Sanchez Almeida Sz Lites (2000). 
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Fig. 9. — Difference between the Stokes I profiles synthesized with and without magnetic fields (J — /non-magnetic)- The difference has 
been modeled as the difference of two Gaussian functions (the dashed line) whose different widths correspond to an apparent magnetic field 
of the order of 220 G. (See text for details.) 
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Fig. 10.— Hanle depolarization produced by the dynamo simulations. The depolarization of Sr i A4607 A is represented versus the height in 
the atmosphere where the depolarizing collisions are evaluated. The two types of line represent different scalings of the dimensionless magnetic 
field strength: the one used along the paper (the solid line), and a factor two smaller (the dashed line). The shaded region corresponds to 
measurements of the depolarization close to the solar disk center (see Faurobert-Scholl et al. 1995). The segment points out the range of 
heights that we ascribe to the dynamo simulation for the synthesis of Fe I 6302.5 A. 
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Fig. 11. — Spatial resolution required to detect a given fraction of the magnetic flux density present in the numerical simulation. The 
diameter of a diffraction limited telescope that yields the required spatial resolution is given in the second axis of ordinates. The solid line 
represents observations without noise whereas the dashed line corresponds to a noise of 20 G, equivalent to a degree of polarization of some 
0.3% . Angular resolutions are in Mm and telescope diameters in m. 




